Improved Lewenstein model for high-order harmonic generation of atoms and 

molecules with scattering wavefunctions 
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We demonstrate a simple method to improve the Lewenstein model for the description of high- 
order harmonic generation (HHG). It is shown that HHG spectra can be expressed as the product 
of a returning electron wave packet and the photo-recombination cross sections, where the former 
can be extracted from the Lewenstein model. By replacing plane waves with scattering waves in 
the calculation of recombination matrix elements, we showed that the resulting HHG spectra agree 
well with those from solving the time-dependent Schrodinger equation. The improved model can be 
used for quantitative calculations of high harmonics generated by molecules. 



43 . 
i 



c3 

1*3 1 

o ■ 

r>3 ■ 
>Y 

43 : 



> 
(N 
0\ 



O 

o 



PACS numbers: 42.65.Ky, 33.80.Rv 

When an atom is subjected to a strong driving laser 
field, one of the nonlinear response processes is the gen- 
eration of high-order harmonics. In the past decade, 
high-order harmonic generation (HHG) has been used 
extensively for the generation of single attosecond pulses 
0) S Si and attosecond pulse trains thus opening 
up new opportunities for attosceond time-resolved spec- 
troscopy. HHG is understood using the three-step model 
[H> @] _ first the electron is released by tunnel ioniza- 
tion; second, it is accelerated by the oscillating electric 
field of the laser and later driven back to the target ion; 
and third, the electron recombines with the ion to emit 
a high energy photon. A semiclassical formulation of the 
three-step model based on the strong field approximation 
(SFA) is given by Lewenstein et al p| . In this SFA model 
(often called Lewenstein model), the liberated continuum 
electron experiences the full effect from the laser field, but 
not from the ion that it has left behind. In spite of this 
limitation, Lewenstein model has been widely used for 
understanding the HHG by atoms and molecules. Since 
the continuum electron needs to come back to revisit the 
parent ion in order to emit radiation, the neglect of the 
electron-ion interaction is rather questionable. Thus, var- 
ious efforts have been made to improve upon the Lewen- 
stein model, by including Coulomb distortion 0, @], or 
by using Ehrenfest theorem 9]. These improvements, 
however, still do not lead to satisfactory agreement with 
exact calculations, nor do they reveal how the target af- 
fects the HHG spectra. 

According to the three-step model, the last step of 
HHG is the recombination process. Taking this model 
seriously, one can ask if the HHG yield can be written 
as the product of a returning electron wave packet with 
the photo-recombination cross section, and if the electron 



wave packet is largely independent of the targets. These 
two assumptions are at the heart of the tomographic pro- 
cedure of Itatani et al. [10]. Their validity has been es- 
tablished for N2 and O2 molecules where the HHG spec- 
tra were calculated using the Lewenstein model and the 
photo-recombination cross sections were calculated using 
the plane- wave approximation (PWA) [ll[ . The indepen- 
dence of the returning wave-packet on the atomic rare gas 
targ ets has recently been investigated by Levesquc et al. 
[12 ] , where the dipole matrix elements (or recombination 
cross sections) were also calculated within the PWA. 

More recently, using HHG spectra obtained from 
numerical solutions of the time-dependent Schrodinger 
equations (TDSE) for atoms in an intense laser field, 
Morishita et al. [l| showed that the "exact" TDSE re- 
sults can be expressed as 



S(u) = W(E) 



da(uj) 
dQi. 



(1) 



where da/dttk is the "exact" photo-recombination (dif- 
ferential) cross section. Here W describes the flux of the 
returning electrons, which we will call a "wave-packet" . 
The electron energy E is related to the emitted photon 
energy u> by the relation E = k 2 /2 = u — I p , with I p be- 
ing the ionization potential of the target (atomic units are 
used throughout this paper, unless otherwise indicated). 
The method for solving the TDSE and the calculation of 
da/dflk are described later. 

In this Letter, we will first show that if the HHG yield 
calculated within the Lewenstein model is written as 



S SFA {u J ) = W SFA (E) d(jPW{uj) 



dfl 



k 



(2) 
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then the electron "wave packet" W SFA will be largely 
identical to W TDSE obtained from Eq. (1) by solving the 
TDSE in the same laser pulse. Here da PW /dClk is the 
recombination cross section calculated the plane wave ap- 
proximation. In the improved Lewenstein model, or the 
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FIG. 1: (Color online) Comparison of the electron "wave 
packets" extracted from the exact HHG spectra of Ne, ob- 
tained by solving the TDSE (solid black line), and from the 
Lewenstein model (dashed blue line) . Also shown is the TDSE 
result for a scaled H (dotted red line) . The laser pulse is with 
duration of 10.3 fs, peak intensity of 2 x 10 14 W/cm 2 and 
mean wavelength of 1064 nm. 



scattering wave based strong-field approximation (SW- 
SFA), we propose that the HHG be calculated by 



S 



SW-SFA 



(w) = W SFA (E) 



da (to) 



(3) 



In Fig. 1, we show the wave-packets extracted from 
the HHG spectra of Ne (solid black line) and from scaled 
H (dotted red line) by using Eq. (1), where the "exact" 
HHG spectra are obtained from the solution of the TDSE 
and the photo-recombination cross sections are calcu- 
lated with "exact" scattering wavefunctions. The effec- 
tive nuclear charge of the companion H is chosen such 
that it has the same Is binding energy as the Ne(2p) 
state. We use a laser pulse with duration (FWHM) of 
10.3 fs, peak intensity of 2 x 10 14 W/cm 2 and mean wave- 
length of 1064 nm. Clearly, the two wave packets agree 
very well. Additional evidences of such agreement can be 
found in our recent work where Eq. (1) was also uti- 
lized to extract very accurate photo-recombination cross 
sections for a variety of atomic targets. The extracted 
electron wave packet W SFA is shown as the dashed blue 
line. This result, which has been normalized to that of 
the TDSE near the cutoff, agrees reasonably well with 
the TDSE result. We have also performed calculations 
for different atoms and with different laser parameters to 
confirm this conclusion. 

Having established the good agreement between W SFA 
and W TDSE , we now present the main results in Fig. 2. 
Here we show the HHG spectra calculated using the di- 
rect solution of the TDSE (solid red lines), Lewenstein 
model (dotted blue lines), together with the SW-SFA 
model (solid black lines) for Ar, Xe, and Ne. For Ar 
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FIG. 2: (Color online) Comparison of the HHG yields ob- 
tained from numerical solution of the TDSE (solid red lines), 
the Lewenstein model (dashed blue lines), and the improved 
Lewenstein model (solid black lines) for Ar (a), Xe (b), and 
Ne (c). 



and Ne, the laser pulse has peak intensity of 2 x 10 
W/cm 2 and mean wavelength of 800 nm. The laser du- 
ration (FWHM) is 10 fs for Ar and 20 fs for Ne. For 
Xe, the corresponding parameters are 5 x 10 13 W/cm 2 , 
1600nm, and 7.8 fs, respectively. The HHG yields for 
Ar are shifted vertically in order to show their detailed 
structures. For Ne and Xe, the SFA and SW-SFA re- 
sults are normalized to the TDSE results near the cutoff, 
i.e., close to 3.2t/ p + I p , where U p is the ponderomotive 
energy. 

First we compare the SFA and the TDSE results. 
Clearly, the Lewenstein model agrees reasonably well 
only in some range near the HHG cutoff for all three 
cases. For Ar, they agree well for a broad range from 
70 eV down to 40 eV. Near 40 eV, the slope of the 
TDSE curve changes and the two results start to devi- 
ate. The discrepancy reaches more than two orders of 
magnitude near the ionization threshold (15.76 eV). For 
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FIG. 3: (Color online) Photo-recombination cross sections of 
Ar (a), Xe (b) and Ne (c), obtained by using exact scatter- 
ing wavefunctions (solid curves) and within the plane-wave 
approximation (dashed curves). 



Xe, the Lewenstein model underestimates the HHG yield 
by up to one or two orders of magnitude, begining just 
below the cutoff. The inaccuracy of the SFA for lower 
harmonics is usually interpreted as due to the neglect of 
the interaction between the target ion and the continuum 
electron in the three-step model. While this interpreta- 
tion is qualitatively correct, it does not explain why in 
some cases the agreement is better. For example, see 
Fig. 2(c) for Ne. Here the HHG yield from the Lewen- 
stein model disagrees with the TDSE results only in the 
range from 28 eV to 45 eV. 

We next examine the prediction of the SW-SFA model. 
From Fig. 2 it is clear that it gives much better agree- 
ment with the TDSE results for all three atoms shown. In 
some cases, the discrepancy has been reduced from two 
orders of magnitude in the SFA to about a factor of two 
in the SW-SFA. Since the recombination cross sections 
calculated using "exact" scattering wavefunctions are in- 
dependent of the electromagnetic gauges used, unlike the 
SFA, the SW-SFA model has no gauge-dependence issues. 
There are some "unphysical" sharp spikes in Fig. 2, near 
22 eV in Ar, 17 eV in Xe, and 32 eV in Ne. They are due 
to the "Cooper minima" where the photo-recombination 
cross sections da PW / dSl k vanish (see Fig. 3). 

In order to understand the nature of the improvement 



in the SW-SFA, we now analyze the photo-recombination 
cross sections. In Figs. 3 (a), (b), and (c) we compare 
the exact results, da/dfl k , with da PW / dfl k for Ar, Xe, 
and Ne, respectively. The figures show that, by approx- 
imating scattering waves by plane waves, the calculated 
dipole matrix elements are quite inadequate for electron 
energies below 100 eV. In particular, in Ne there is a spu- 
rious "Cooper minimum" near 32 eV if the plane wave 
is used. This minimum is totally absent for the whole 
energy range in the exact result. The Cooper minima 
are known to be related to the zeros in the dipole ma- 
trix elements and they are absent for Ne since the 2p 
wavefunction is nodeless We note, from Fig. 3(b), 

the position of the Cooper minimum for Ar predicted by 
using the PWA is shifted lower by about 20 eV as com- 
pared to the exact result. For the Xe case, the shift is 
even larger (the "exact" Cooper minimum, which is near 
80 eV, is not clearly seen in Fig. 3(b)). 

The theoretical method for the numerical solution of 
the TDSE has been described previously 15, Il6j |. The 



electric field of the laser pulse is written in the form 
E(t) = Eoa(t) cos(uit + <f>), with the envelope given by 
a(t) = cos 2 (7rt/r), where r is 2.75 times the FWHM of 
the laser pulse. All the atoms are approximated by a 
one-electron model potential V(r) which has the form of 
V(r) = V s (r) — 1/r. The parameters in the model po- 
tential are chosen such that the binding energies of the 
ground state and the first few excited states are close to 
the experimental values. We then propagate the initial 
ground state in time using the split operator method. 
The laser parameters were chosen such that the deple- 
tion effect is small. Typically the ground state survival 
probability is more than about 95%. 

Photo-recombination cross section is determined as 



da (to) 
dflh 



27rfc 



*l*o)| 



(4) 



where \^ + ) and |<E>o) are the scattering wavefunction and 
the ground state wavefuntion, respectively, and £l k = 
is the direction of fe. The returning electron 
that contributes most to the HHG is the one that prop- 
agates along the laser polarization direction. Therefore, 
the relevant cross section is for k that is parallel to z-axis 
(the laser polarization axis), that is, 8 k = or tt. 

For each model potential V(r), continuum scatter- 
ing wavefunctions are calculated numerically, using the 
partial-wave expansion 



i=o 



k 



i 



Y{ m ie k ,(j> k )Y lm {e,cj>) (5) 



where 8i(k) is the partial-wave phase shift. With this 
scattering wave function, "exact" photo-recombination 
cross section is then obtained by using Eq. 4. Note that 
for atoms which have po ground states, only / = 0, 2 and 
m = will contribute to the dipole matrix element. It 
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molecules, we expect that accurate HHG spectra from 
molecular targets can be calculated. For molecules, ac- 
curate solutions of the TDSE are very time-consuming, 
if not impossible. Using the Lewenstein model, the 
HHG as well as the dipole matrix elements can be eas- 
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FIG. 4: (Color online) The electron wave-packets, extracted 
from the HHG spectra for Ne and Ar under the same laser 
field. The inset shows the wave-packets in the linear scale 
with the Ne data shifted horizontally by -1.2 eV. The HHG 
are obtained within the Lewenstein model with 1064 nm laser, 
peak intensity of 2 x 10 14 W/cm 2 , duration of 50 fs. 



is interesting to note that for all the systems considered 
here, the shape of da/dilk as functions of photon energy 
are quite close to the total (integrated) cross sections, 
obtained by averaging over m = —1,0,1 initial states. 
In the PWA, a plane-wave (2tt)~ 3 ^ 2 exp(ifcr) is used to 
obtain d<r PW / d£l}~. 

So far we have extracted photo-recombination cross 
sections from Ar, Xe, and Ne by using companion atoms 
with nearly identical ionization potentials, as compared 
to the atoms under consideration This is used to 

greatly simplify the analysis. In the following, we will 
show that this requirement can be relaxed. As the wave 
packets obtained within the Lewenstein model agree rea- 
sonably well with the TDSE results, for our purpose we 
only use the Lewenstein model in the subsequent analy- 
sis. In Fig. 4, we show the comparison of the wave packet 
from Ne (I p = 21.56 eV) and Ar (I p = 15.76 eV) in the 
same 1064 nm laser pulse with peak intensity of 2 x 10 14 
W/cm 2 , duration (FWHM) of 50 fs. Clearly, the two 
wave packets lie nicely within a common envelope. This 
can be seen even more clearly in the inset, where the 
data are plotted in linear scale with the Ne data shifted 
horizontally by -1.2 eV. We note that the small details 
below 33 eV also agree well. This conclusion is also con- 
firmed by our calculations using the TDSE with different 
laser parameters and with other atoms. This supports 
that the independence of the wave packet on the target 
structure can also be extended to systems with different 
ionization potentials. 

In this Letter, we have shown that the improved 
Lewenstein model is capable of describing HHG spec- 
tra close to the accuracy that can only be achieved by 
solving the TDSE directly. Extending this model to 



ily calculated [17j, [18( to extract the wave packet. Sim- 
ilarly, theoretical tools for the calculation of molecular 
photo-recombination (or photoionization) cross sections 
are quite mature, and general computer codes are avail- 
able even for large molecules (Tij]. Thus it is expected 
that HHG spectra generated by molecules can be cal- 
culated using the SW-SFA model readily. Finally, we 
mention that the present analysis has not yet addressed 
the phases of the harmonics. Further studies would be 
needed to disentangle the phases due to the laser and the 
target separately. This step is needed in order to incorpo- 
rate propagation effect of the medium on the generated 
high harmonics. 
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